Efficient multiparticle entanglement via asymmetric Rydberg blockade 
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We present an efficient method for producing TV particle entangled states using Rydberg blockade 
interactions. Optical excitation of Rydberg states that interact weakly, yet have a strong coupling to 
a second control state is used to achieve state dependent qubit rotations in small ensembles. On the 
basis of quantitative calculations we predict that an N = 8 Schrodinger cat state can be produced 
with a fidelity of 84% in cold Rb atoms. 

PACS numbers: 03.67.Bg, 32.80.Qk, 32.80.Ee 



Entanglement lies at the heart of quantum informa- 
tion processing and is also a valuable resource for extend- 
ing precision measurements beyond bounds set by clas- 
sical statistics. Recent years have seen a steady progres- 
sion towards entanglement of larger and larger objects. 
Although macroscopic ensembles have been successfully 
cntangledpj, the entanglement achieved per atom was 
very low. Maximally entangled cat states of six atoms, 
as well as "W" states of eight atoms have been achieved 
in groundbreaking experiments with cold ions0, In 
this letter we introduce an efficient technique for gener- 
ating maximally entangled states which is applicable to 
any system which supports asymmetric state dependent 
blockade interactions. We give quantitative estimates 
for the preparation fidelity for entanglement of the clock 
states of Rb atoms using Rydberg blockade, which may 
enable improvement in the accuracy of an atomic clock. 




FIG. 1: (color online) Level scheme (left) and sequence of 
operations for cat state generation (right), is the effective 
Rabi frequency coupling states |0), |1). 



Consider the situation shown in Fig. 1 where N atomic 
quits, each with basis states |0), |1), are confined in a 
volume V. We assume states |0), |1) are weakly interact- 
ing over time scales of interest but can be transferred 
to additional interacting states |s),|p). Single particle 
excitations of \s) are allowed but there is a large en- 
ergy gap U ss — hA ss which blocks two-particle excita- 
tions. States |s), \p) are also strongly interacting with a 
large gap U sp = hA sp , however states \p) interact weakly 
with each other so that the two-particle interaction en- 
ergy U pp = hApp satisfies A pp < A sp , A ss . 



With the above resources N atom entangled states can 
be synthesized in a few interaction steps by the following 
protocol. We first prepare the N atom product state 
\tjj) = |0, 0, ...0). The ground state |0) is coupled to \s) 
with an interaction Hamiltonian TL\ such that the Rabi 
frequency (from now on we put h = 1), given by Q s /2 = 
(s|Wi|0), satisfies |fi s | -C A ss . In step (i) we apply 7i\ 
to all atoms for a time t\ = 7r/(2\/]V|f2 s |) to create the 



entangled state 




A' 



^|0,0,s«...0) + |0,0,...0) 



(1) 



3=1 



We then invoke a second interaction Hamiltonian 7i 2 = 
+^2i with corresponding Rabi frequencies f2 p o/2 = 
(p|^2o|0), fipi/2 = (p|?^2i 1 1) , and the same detuning 
Ao on both transitions, see Fig. 1. For simplicity we 
will assume f2 p o = f2 p i = f2 p = |O p |. After a definite 
interaction time of t 2 = \f2it/Vlp in the resonant (Ao = 
0) case, and t% = 27rA /fi p in the non-resonant (A ^ 
ftp) case, H2 induces a transfer from |0) to |1) in all the 
atoms via the Rydberg state \p), unless this process is 
blocked by population in the Rydberg |s) state. In the 
limit where A pp <C £1 <S A sp step (ii) transforms (TTJ) into 
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We finish in step (hi) by applying —Tii for a time 2ti to 
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reverse the first excitation step, giving 



1 



V2 



(|0,0,...0) 



1,1, 



•1» 



(3) 



which is a N atom maximally entangled state. We see 
that, independent of N, only three preparation steps 
are needed. We show below that the requirement of 
strong and state dependent asymmetric couplings may 
be satisfied by dipole-dipole interactions of Rydberg 
atoms 0,(1]. Rydberg blockade effects have now been ob- 
served in a number of experiments in both many-body 
and single atomQ settings. The Rydberg blockade has 
been suggested previously as a route to multi-particle 
entanglement HQ and recent work from Midler et al. [Toj 
is based on ideas closely related to those presented here. 
There are, however, significant differences including our 
use of blockade in step (i) of the above protocol which re- 
moves the need for separately addressing a control atom. 
This allows all atoms to reside in one ensemble which 
provides a better geometrical scaling of the interactions. 

The fidelity with which state ([3]) can be prepared in 
an experiment depends on the strength of the blockade 
interactions and the degree to which the couplings are 
asymmetric. We proceed by estimating the effective Ry- 
dberg interaction strengths A sp ,A P p. The interaction 
between atoms in the "control" state \s) and the "target" 
state \p) is of resonant dipole nature between two-atom 
states \sp), \ps). The interaction between classical dipoles 
is anisotropic and has a zero when the angle between the 
dipoles is 9 = cos _1 (l/v / 3). This anisotropy would lead 
to unacceptable errors in the present setting. However, 
for small external fields the atomic Zeeman states are 
degenerate and the interaction couples states with dif- 
ferent m si m p quantum numbers. This leads to a finite 
interaction strength A sp at all angles. It can be shown 
that the behavior corresponding to interaction of classi- 
cal dipoles with angular zeroes is recovered by applying 
a large magnetic field that selects a single pair of Zee- 
man states. The resonant interaction has a 1/R 3 scaling 
so we can write A sp (i?) = A sp (d)(d/ R) 3 where d is a 
characteristic length scale that we will set equal to the 
smallest interatomic separation d. The target-target in- 



teraction A pp is due to a Forster process |12| . We will 



choose states and values of d such that we are working 
in the van der Waals limit of this interaction which gives 
the distance scaling A pp (i?) = A pp (d)(rf/i?) 6 . The en- 
ergy shift of each atom is thus dominated by its nearest 
neighbors in the ensemble. 

bmce A sp ~ 1/i? 3 and A pp ~ 1/i? 6 the condition of 
strongly asymmetric Rydberg interactions can be readily 
met by choosing R sufficiently large. The asymmetry is 
maximized for small n since the resonant dipole allowed 

4 while the 



sp 



interaction between | s) , \p) scales as A 
second order Forster process leading to A pp scales as 
. The lower limit on n is set by the blackbody 



A 



,11 



pp 



limited spontaneous emission lifetime r p ~ n . 

We have searched for parameters in the small n regime 
that minimize the error in creation of state ^ by per- 
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FIG. 2: (color online) Calculated blockade shifts as a function 
of the angle of the molecular axis for 87 Rb, \s) = |41s 1 / 2 ,m = 
1/2), \p) = |40p 3/2 ,m = 1/2) with magnetic field B = 10~ 7 
at d = 3 fim, 

the indicated Forster channels. 
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forming extensive numerical studies of the interaction be- 
tween different Rydberg states in Rb drawing on the ex- 
position of the dominant Forster channels given in [l3|. 
As shown in Fig. [2] we find large interaction asymmetries 
using n s — 41 s 1 / 2 states for \s) and n p = 40 P3/2 states 
for \p) at a lattice spacing of d = 3. /im. In calculating 
the \p) state interactions A pp we have included all chan- 
nels of the form npnp — ► n'sn" s, n'sn"d, n' dn" d. We see 
from the figure that A sp /A pp > 150 for all angles . In 
addition to multiparticle entanglement this large asym- 
metry will also facilitate implementation of a three-bit 
Toffoli gate[il]|. 

We consider an implementation with a cubic lattice 
of spacing d occupied by one atom per site inside a 
sphere of radius Rq. A protocol for preparing a lattice 
with this type of spatially localized occupation was de- 
scribed by us recently in Ref. [l5j . The angle depen- 
dent peaks in A pp are not of particular concern since the 
cubic lattice can be oriented to avoid the corresponding 
angles. For simplicity we have characterized the angle av- 
eraged interactions by A sp (d) = j^ 2 d9 A sp (d, 9) sin(0) 
and similarly for A pp . Performing the integrations we 
find (A sp , A pp )/27T = (14.4,0.019) MHz. 

In order for step (i) of the entanglement protocol to 
be effective it is also necessary that the \n s = 41,Sx/ 2 ) 
states exert a strong blockade over the entire sphere of 
radius i?o. At d = 3 /mi we find A ss (d) = 3.7 MHz and 
the interaction is essentially isotropic T3]. This is larger 
than A pp but the interaction strength decreases as 1/i? 6 
and is insufficient for strong blockade over a sphere with 
Ro > d. We note that this difficulty can be readily solved 
as follows. The first Rydberg excitation step (i) is made 
to a level \s') that has a large value of n and provides 
strong blockade over the entire ensemble. The level \s') 
is then transferred to \s) using a two-photon transition 
which prepares the state of Eq. (P) even though A ss 
may be small. The additional transfer steps are then 
run backwards in step (iii) . The possibility of performing 
these additional steps allows us to ignore the small errors 
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associated with blockade of the control state. 

To examine the fidelity of our scheme it is instructive 
to recall how the errors scale for a two-atom Rydberg 
blockade phase gate. A phase gate between a control 
atom (c) and target atom (t) involves the stepsQ: i) it 
pulse |l) c -> i\r) c , ii) 2ir pulse |l) t -> i\r) t -> -|l)t, 
iii) 7r pulse i|r) c — > — |l) c . Assuming ground to Ryd- 
berg state oscillations can be driven with high accuracy, 
as has been demonstrated in recent experiments[ll|, the 
dominant errors come from imperfect blockade in step 
ii) with error E^\ ~ f2 2 /A 2 (A is the dipolc-dipole in- 
teraction shift) and spontaneous emission of the control 
atom with error E sc ~ 1/fir where r is the Rydberg 
state spontaneous lifetime. The sum of the two errors is 
minimized for Q ~ A 2 / 3 /r 1 / 3 which leads to a gate error 
that scales as0 E = E u + E sc ~ l/(Ar) 2 / 3 . Typical 
numbers for atoms separated by R ~ 5 /im and Rydberg 
principal quantum number n ~ 75 are A ~ 27r x 5 MHz 
and r ~ 200 ^s leading to errors E < .01. 

Similar error estimates apply to the above entangle- 
ment protocol, but in addition to the imperfect blockade 
and the atomic spontaneous decay, we must also take into 
account the undesired interaction between atoms in the 
Rydberg \p) transfer states. The variation in this inter- 
action comes, on the one hand, from considering atoms 
with all nearest neighbors present relative to atoms at 
the edge of the ensemble with fewer neighbors and, on 
the other hand, from the quantum mechanical spreading 
of the occupancies of the Rydberg state \p) in the time 
evolving many-atom superposition states. 

Returning to Fig. [TJ we expect process ii) to be fast 
and most sensitive to the blocking interacting A sp , but 
also most sensitive to the A pp shifts, when Ao = 0. 
We will treat the spontaneous decay and the imper- 
fect blocking due to finite A sp as independent errors 
on each atom. The spontaneous emission error during 
the transfer between |0) and |1) via \p) is readily deter- 
mined from the average population of the Rydberg state 
to be E sc = ^pf-. The states |0) and |1), coherently 
coupled to \p) with equal Rabi frequencies Cl p , can be al- 
ternatively treated in the basis of the uncoupled, "dark" 
state \d) = (|0) - |l))/\/2 and the coupled, "bright" state 
1 6) = (|0> + |l))/\/2 with Rabi frequency V2Q p . A 2tt 
rotation on the \b) — \p) transition yields a sign change 
on the bright state, equivalent to the desired net it rota- 
tion between the |0) and |1) states. This motivates the 
definition of f2 = as an effective Rabi frequency 

of the oscillation between |0) and |1), and in the limit 
of O l/r p the populations of states |0), |1) after the 
transfer pulse (ii) are Pq = I — P±, 

Pi = [n 12 - n 2 + n 2 cos(7rO'/o) 

- 0' 2 cos(7rA/2^)cos(7rrj72r2) 

- AO'sin(7rA/2Q)sin(7rfi72fi)] /(2fl' 2 ), (4) 

il' = \/A£l 2 + A 2 , and A is the interaction induced de- 
tuning. The error due to imperfect blocking by the con- 
trol state \s) is found by inserting A — > A sp 3> ft and this 



leads to Eu ~ ^- \ > where we assume an average 
over atom pairs in the ensemble. 

Equation (TJJ also allows an estimate of the transfer 

error on each atom scaling as E tI cx A ^ 2 d due to the 
A pp interaction terms, when the transfer is not blocked, 
but we will evaluate this error taking into account the 
full many-atom correlations in the quantum state. Our 
task is to determine the effect of the interaction: 

V='£AV p (\p i )(p i \)®(\p j )(p j \), (5) 

i,j>i 

where \pi){pi\ is the projection operator of the i fh atom 
on the Rydberg state \p), and A l J p is the interaction 
energy for a given (ij) pair of Rydberg excited atoms, 
depending on their spatial separation. We will deter- 
mine the effect of this interaction by first order per- 
turbation theory, in the interaction picture with respect 
to the ideal gate operation due to the Hamiltonian TI2 
transferring the atoms between state |0) and |1). That 
Hamiltonian is readily diagonalized for each atom, 7i 2 z = 
X)m w ml TO j)( TO ili where both the energies oj m and the 
states \rrii) of the i th atom are analytically known (and 
the same for all atoms) , and the corresponding ideal time 
evolution operator U^it) of the entire atomic ensemble is 
thus also known. The time evolution in the interaction 
picture, due to the Rydberg interaction is given to first 
order by the expression: 

U T (t) = I-iJ2 A ppYl C ™* C ™3 C *m> C *m> 
i,j>i mi,m,j m^m'j 

t(u m . -\-uJ m . — u / — u? / )t 

e 1 1 i i — 1 

x(\m i )(m' i \)®(\m j ){m' j \)® (g) I k , (6) 

where c mi = (rrii\pi) comes from the expansion of the 
Rydberg states on the eigenstates of Jin. 

The perturbation leads to an erroneous change of the 
state in the interaction picture, and one readily observes 
that the squared norm of the erroneous component has 

A 2 

the expected -ffi scaling. We have evaluated © by sum- 
ming over all atom pairs and assuming the initial state 
with all atoms in state |0), and we find for different atom 

numbers the transfer error Btr,iv = «at • 02 ■ We have 
carried out calculations with a pair of atoms, separated 
by the distance d, four atoms located at the corners of a 
square with side length d, and 8 atoms located at the cor- 
ners of a cube with side length d. With the R~ 6 scaling of 
the interaction with distance, we find (a^\oi£\c^p) = 
(0.299,0.72,9.39), while for a distance independent cou- 
pling, (4 0) ,4 0) ,4 0) ) = (0.299,3.82,36.8). The accu- 
racy of the perturbation theory results has been verified 
by direct numerical solution of the N atom Schrodinger 
equation for N < 6 . 
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FIG. 3: (color online) Error of the N — 8 cat state calculated 
from Eq. using the parameters of Fig. 2, and r p = 57 /us 
for n p = 40. 



to the state of Eq. @ gives in this case a spontaneous 
emission error per atom of E sc = -| A * r . In the limit 

of D, <C Ao the population of the states |0), |1) after the 
transfer pulse are Pq = cos (nta/2), Pi = sin 2 (m 2 /2) 
with Q = Qp/2(A + A sp ) when the \s) state is supposed 
to block the transition. The blockade error in the tar- 
get state probability for each atom is then determined 

as Em = x (xf~) ■ When the transition is not blocked, 
we shall use our perturbative expression ([5]), which in 
the non-resonant case yields an expression of the form 

Ett n = Pn ■ — ^3—- With the same arrangement of 4 
and 8 atoms as above we find with the R 

?(6) o(6) 



interaction, 



(fll ,(3 S ) = (15.6,113), while for a distance indepen- 
dent coupling, (/?4°\/3g^) = (53.7,308). In this case, the 
errors add to 



The approximate cubic growth of the transfer error 
with the number of atoms qualitatively agrees with an 
error amplitude on each atom scaling linearly with the 
number of perturbing atoms. For larger ensembles, dis- 
tant neighbors do not contribute to the error, and we 
expect a transition to a linear dependence with N. 

Adding together our error contributions, we find the 
total error on the N-atom state, 



E N = N 



7T 1 
4 ilTn 



7T 

T 



r> 2 

A 2 



KM) 



(7) 



o 2 



As can be seen in Fig. [3] the contributions to the error 
depend in different ways on the Rabi frequency. We find 
for the case of N = 8 atoms a minimum error of E$ = 
0.16 at S1/2tt = 0.30 MHz. The 8 atom cat state can thus 
be prepared with a fidelity of ~ 0.84. 

We have also performed the calculation assuming 
Ao 3> App(d). The it pulse in transfer step (ii) leading 



E N = N 



1 



2 A r p 



7T 

T 



A 2 
A 2 



KM) 
K ' 



(8) 



As the (3m coefficients are substantially larger than the 
corresponding a at the non-resonant transfer case has a 
lower fidelity than for the resonant case. 

In summary we have presented a technique for prepar- 
ing multi-atom maximally entangled states using a three 
step sequence. A detailed analysis of asymmetric Ry- 
dberg interactions in Rb atoms shows that 8 atom cat 
states can be prepared with reasonably high fidelity. Sim- 
ilar results, not presented here, have been found for the 
case of Cs. Straightforward modifications to these ideas 
can be used for two-atom CNOT gates that do not re- 
quire single qubit rotations [17|, and in other physical 
settings where blockade interactions are available such 
as Coulomb or Pauli blockade of quantum dotsflSfl. or 
molecular interactions with superconducting qubits[19(. 

This research was supported by the NSF, ARO, and 
IARPA. 



[1] B Julsgaard, A. Kozhekin, and E. S. Polzik Nature 413, 
400 (2001). 

[2] D. Leibfried, et al., Nature 438, 639 (2005). 

[3] H. Haffner, et al., Nature 438, 643 (2005). 

[4] D. Jaksch, et al., Phys. Rev. Lett. 85, 2208 (2000). 

[5] M. D. Lukin, et al. Phys. Rev. Lett. 87, 037901 (2001). 

[6] D. Tong, et al., Phys. Rev. Lett. 93, 063001 (2004); K. 

Singer, et al., Phys. Rev. Lett. 93, 163001 (2004); T. 

C. Liebisch, et al., Phys. Rev. Lett. 95, 253002 (2005); 

Erratum, ibid 98, 109903 (2007); T. Vogt, et al., Phys. 

Rev. Lett. 97, 083003 (2006); R. Heidemann, et al., Phys. 

Rev. Lett. 99, 163601 (20 07). 
[7] E. Urban, et al., larXiv:0 805.0758 (to appear Nature 

Physics); A. Gaetan. et al.. larXiv:0810.2960l 
[8] R. G. Unanyan and M. Fleischhauer, Phys. Rev. A 66, 

032109 (2002). 

[9] D. M0ller, L. B. Madsen, and K. M0lmer, Phys. Rev. 
Lett. 100, 170504 (2008). 
[10] M. Miiller, I. Lesanovsky, H. Weimer, H. P. Biichler, and 



P. Zoller, larXiv:0811.lT55l 
[11] T. A. Johnson, et al., Phys. Rev. Lett. 100, 113003 
(2008). 

[12] T. Forster, Ann. d. Physik 437, 55 (1948). 
[13] T. G. Walker and M. Saffman, Phys. Rev. A 77, 032723 
(2008). 

[14] E. Brion, A. S. Mouritzen, and K. M0lmer, Phys. Rev. 

A 76, 022334 (2007). 
[15] M. Saffman and K. M0lmer, Phys. Rev. A 78, 012336 

(2008). 

[16] M. Saffman and T. G. Walker, Phys. Rev. A 72, 022347 
(2005). 

[17] N. Ohlsson, R. Krishna Mohan, and S. Kroll, Opt. Com- 

mun. 201, 71 (2002). 
[18] M. Friesen, A. Biswas, X. Hu, and D. Lidar, Phys. Rev. 

Lett. 98, 230503 (2007). 
[19] K. Tordrup and K. M0lmer, Phys. Rev. A 77, 020301(R) 

(2008). 



